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^ ■ Abstract 

The recently observed X(3872) resonance, which is difficult to be assigned a conventional cc 
. charmonium state in the quark model, may be interpreted as a molecular state. Such a molecular 

^sO ' state is a hidden flavor four quark state because of its charmonium-like quantum numbers. The 

s-channel one gluon exchange is an interaction which only acts in the hidden flavor multi-quark 
system. In this paper, we will study the X(3872) and other similiar hidden flavor molecular states 
in a quark model by taking into account of the s-channel one gluon exchange interaction. 
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I. INTRODUCTION 



Recently, many new resonances are discoveried experimentally. Many of them have the 
proper quantum numbers of the qq meson states. However, their mass values do not fit the 
conventional qq states in the quark model. Among them, X(3872) was first observed in the 
J/i/jir + ir~ channel by Belle collaboration in 2003 [1], and has been confirmed by CDF 2], 
DO [2I and Babar collaborations 0]. Its quantum numbers are probabely J PC = 1 ++ . The 
corresponding charmonium candidates in the quark model are 2 3 P!(3990) and 3 3 Pi(4290) 
which are 50 ~ 200 MeV above M x = 3872 MeV. 

Many people suggested that X(3872) is mainly a DD* molecular state 
However, to bind the quarks and anti-quarks together in such a four quark state or other 
multi-quark states, we need to introduce new interaction into the quark model. Swanson 
proposed that the X(3872) is mainly a DD* molecule bound by the meson- meson interaction 
derived from the one pion exchange and the quark exchange ?]]. In Wong's work [8], the 
meson-meson interaction is derived from a QED-type effective interaction in terms of effective 



J~T 
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15], further investigations 



charges for quarks and antiquarks. In refs. [lOl . 
based on the molecule assumption were carried out. 

Since the strength of the one pion exchange interaction seems not strong enough to 
bind the DD* molecular state, other authors argued that X(3872) may be a dominant cc 



charmonium with some admixture of DD* 
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17j, ll8|. In ref. [3], after taking into account 



of the sigma meson exchange potential, the interpretation of X(3872) as a loosely bound 
molecular state was further disfavored. 

We should notice that the color structure of a multi-quark state is much richer than that 
of a qq conventional meson state. Unlike the conventional mesons or baryons, the qq and qq 
pairs in a multi-quark state can be in the color 8 C and 6 C representations respectively. Some 
color interactions which have no effects in the qq or g 3 colorless system, may have significant 
contribution in a multi-quark system. So the complete interactions in the quark model can 
be quite different after we take into account of these multi-quark states. 

The s-channel one gluon exchange interaction is an interaction between quark and anti- 
quark of the same flavor which annihilate into a virtual gluon. It has no effect on the conven- 
tional qq mesons but acts in the hidden flavor multi-quark system like the charmonium-like 
moelcular states. In this work, we will investigate the hidden flavor molecular states by con- 
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sidering the s-channel gluon exchange interaction in the quark model. In the next section, 
we will model the potential of the s-channel one gluon exchange interaction starting from 
the non-relativistic reduction. Then the J PC quantum numbers of the molecular states are 
selected by an analysis of the spin dependence of the interaction strength. In sec. IHIl we 
will carry out the numerical calculation of X(3872) and some other charmonium-like states 
as the molecular states. Also we will make prediction about the similar bottomium-like 
molecular states. Finally, we will give a brief summary. 



II. THE POTENTIAL OF S CHANNEL ONE GLUON EXCHANGE INTERAC- 
TION 



In our work, we will use the Bhaduri quark model which is a rather simple non-relativistic 
quark potential model. In the Bhaduri model, the hamiltonian can be written as 20J 



» = I> + S-lE( f <- F ' V S + • FA ■ SjVf) . 



(i) 



i<j 



rrii are constituent quark masses and Ff = ~A? (c = 1,...,8) are the well-known SU C (3) 
Gell-Mann matrices. Apart from a constant, here the central potential is the usual one 
gluon exchange coulomb potential plus the linear confinment: 



Vg — ■ -f-.U () . 



(2) 



|r*j — rA is the distance between quark i and j. The color- magnetic interaction reads 



V; 



ss 



4k 



miiTLj rfirij 



-rtj/ro 



where the 5-interaction has been smeared smoothly with the prescription 

1 _-../»... 



(3) 



(4) 



The model parameter values are 



k = 102.67MeVfm, 
M = 913.5MeV, 
fn u = f^d = 337MeV, 
m c = 1870MeV, 



a = 0.0326(MeV~ 1 fm)i 

r = 0.4545fm 
m s = 600MeV, 
m b = 5259MeV. 
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FIG. 1: The s-channel one gluon exchange 




p', s> 



P, s -p, s 



The s-channel one gluon exchange interaction (SOGE) takes place when a qq pair anni- 
hilates into a virtual gluon (Fig. CE|). The non-relativistic reduction of the potential is: 



(5) 



where G(4m 2 ) is the one gluon exchange amplitude. The last 5(fc, fj) represents that the 
quark and the anti-quark must have the same flavor. Clearly, the two factors in the brackets 
mean that this interaction only occurs when the qq pair is in the color octet and with 
spin S = 1. In our model calculation, the 5- interaction should be smeared smoothly with 
the same prescription used in the color-magnetic interaction in the same quark model. 
However, the QCD amplitude G(4m 2 ) is taken value in the timelike region, where the QCD 
behavior is still not clear at present time. In the spacelike region, this amplitude is well 
known from the gluon propagator in perturbative QCD and it reads 

GUM 2 ) = ^ for q 2 < 0. (6) 

However, in our study of the hidden flavor molecular states, in order to provide an attrative 
interaction to favor the formation of the bound states, we need that 

G(4m 2 ) > 0, 

which means the above formula from perturbative QCD should not be directly used in the 
timelike region. This change of sign is first suggested in ref. [21] in the study of irir and Kir 



scattering. Following ref. 



21], we assume that 



G(4m 2 ) = -fGUMm 2 ) = /— f , (7) 



where / is an introduced strength factor. In our model, after the ^-function being smeared 
smoothly, the SOGE potential turns to be: 

vr * = -{(l + |^.^) (| + ^. 8,) ^e-^ (/ „ f]) . (g) 

First we observe that the potential y SOGE is proportional to m^ 2 , so the interaction 
mainly takes place between the light qq pairs of q — u, d, s in any multi-quark system. 

Next, we will analyse the spin dependence of ^/ SOGE in each J PC channel of molecular 
state. For the molecular states we concern, let us assume that their flavor structure is QqqQ, 
where Q = c,b and q — u,d, s. Since we can neglect the K SOGE interaction between QQ, 
the interaction acts only when the qq pair is in color octet and S = 1 as we have mentioned 
before. More specifically, the interaction will be switched on if the spin coupling of the four 
quarks is 

[(QQ)ji(qq)j 2 =i\j- 

The color structure of QQ and qq should be octet obviously. However, the color (re) coupling 
is irrelevant to our analysis here and will not be explicitly presented. Furthermore, we will 
assume no spatial excitation of any quark. Hence all quarks have zero orbital angular 
momentum. Then the quantum numbers J PC can be determined easily for J\ = 0,1. We 
have following four \/ SOGE interaction channels: 

Ji = 0:l + ~; Ji = 1 : ++ ,l ++ ,2 ++ . 

Finally, we can recover the molecular states by the angular momentum recoupling. We 
obtain 



0++ 



i++ 



1+- 



[(QQ)i(qq)i]o = -^-(QqMqQ)o ~ \[(Qq)i(qQ)i]o, 



[(QQ)i(qq)i]i = -^=(Qq)i(qQ)o + ^=(Qq)o(qQ)i, 



(QQ)o(qq)i = \(Qq)i(qQ)o + \(Qq)o(qQ)i + ^=[(Qq)i(qQ)i]i, 
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TABLE I: The list of spin factor of y SOGE interaction strength in a molecular state J PC made of 



two meson states J^ 1 and J 2 2 ■ 


c.c. = charge conjug 


ation. 


J [up 


jPC 


Spin factor 




0++ 


3 
4 




0++ 


1 

4 


i~o- +c.c. 


1++ 


1 


l - _ -c.c. 


1+- 


1 

2 




1+- 


1 

2 


l-i- 


2++ 


1 



• 2++ 

[(W)l(??)l]2=[(Q9)lW)l]2 



Then the factor of interaction strength can be read from the coefficients, which we present 
in table [H We see that the \/ SOGE interaction favors to form the molecular states with 
J PC = l++,2++. 



III. NUMERICAL CALCULATION 

To calculate the molecular states, we use the Rayleigh-Ritz variation principle. The test 
wave function will be taken to be a series of Gaussian basis functions. The Gaussian basis 
functions are often utilized in variational calculations of atomic and molecular problems. 
Recently the method has been also used in the few body system in nuclear and particle 



physics 
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I24j. 

In our case of the QqqQ molecular state, the test wave function of a molecular state 
between two clusters of qq meson states is a series 

^i234(ri2, r 34 , ri2 34 ) = ^ "1234^12 (^12)^34 (r 34 ) exp(-(3[ 2U rl 2U ), (9) 

i 

where 7*1, r 2 , T3 and r 4 are the coordinates of Q, q, q and Q, respectively, = r { — r,-. 
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r i234 is the distance between the two meson clusters 

m Q v x + m q r 2 m q r 3 + m Q r± 

?"i234 = ; ; • (,-WJ 

rriQ + m q m q + m Q 

i J ij( r ij) 1S the meson wave function which is also taken to be a Gaussian function series 

M-rti) = 4 exp(-4-rJ). (11) 

k 

The wave function of a molecular state is determined by the variation principle in two 
steps. We first determine the wave function (1111) of each meson cluster. Then the meson 
cluster functions ipij are fixed in to obtain the wave function of the molecular state and 
their masses. 

To reduce the amount of computation, the parameters f3 l and a 1 in a Gaussian function 
series are determined in two steps by one-dimensional minimization. We first determined a 
average (3 value using a single Guassian function. Then a set of 2N + 1 elements is 



generated from scaling the j3 value up and down by a scale factor s 

P = Ps'- N (12) 

where i = 0, 1,...,2N. The coefficients a 1 are determined by diagonalization the model 
Hamiltonian in the 2N + 1-dimensional space spanned by these 2N + 1 different Gaussian 
functions. The final values of f3 % and a % and the mass of tetra quark states are determined 
by searching the scale factor s for a minimum of system energy. 

In this way, we have calculated the possible ++ (the combination of two 0~ mesons with 
the spin factor |), 1 ++ and 2 ++ molecular states in table [H We have calculated both the 
charmonium-like states cqqc which will be compared to recent experimental results and the 
bottomium-like states bqqb which can be investigated in further experiments. 

The results of relevant heavy quark Qq mesons are listed in table [III We see that the 
mass values calculated from the Bhaduri quark model deviate at most 30 MeV away from 
the experimental data. The difference between the calculation from the variation method 
with Gussian function series (cal. II) and the exact numerical calculation (cal. I) is less than 
0.5 MeV. So the variation method is an impresssive good approximation for the numerical 
calculation of the Qq conventional mesons. 

The results of the molecular states are given in table [TTT1 and table HVl Here the molecular 
state is presented by its binding energy: 

E b = M 1 + M 2 -M x , (13) 



TABLE II: Mass of Qq mesons. The experimental values are taken from PDG 25|. In calculation 
I, the mass is obtained by solving the Schrddinger equation. In calcalation II, the mass is obtained 
by the variation method with Gaussian test functions using N = 3. 



meson state 


exp. (MeV) 


c al. I (MeV) 


c al. II (MeV) 


£>± 


1869.62 


1885.56 


1886.14 


D° 


1864.84 






£>*(2007)° 


2006.97 


2019.96 


2020.07 


L>*(2010) ± 


2010.27 






Df 


1968.49 


1995.78 


1996.36 


D*± 


2112.3 


2101.2 


2101.4 


S± 


5279.15 


5300.84 


5301.18 


B° 


5279.53 






B* 


5325.1 


5350.3 


5350.5 


B° 


5366.3 


5371.9 


5372.4 


B* s 


5412.8 


5413.3 


5413.6 



and the rms radius (r 2 ) 1//2 , where Mi, M 2 are masses of the two compound mesons, Mx is 
the mass of the molecular state. 

We can see from tables IHII and IIVI that indeed the 1 ++ and 2 ++ molecular states are 
bound preferable since the ^/ SOGE are strong in these two channels. The binding of 1 ++ 
states are slightly deeper and tighter than that of 2 ++ 's. So from our model calculation, the 
1 ++ molecular states should be easily observed in experiments. 

First, let us look at the DWDM sector. If the X(3872) is a molecular state of D and D* , 



the binding energy can be estimated from experimental data|25| 

E b = M D + M D * -M x = 0.3 MeV, 

which is very small. Accounting the uncertainty of approximation in our model calculation, 
we think the reasonable range of / value should be around —0.8 ~ —1.0. If it is small 
/ ~ —0.8, then it may be the only molecular states in the £)(*)£)(*) sector. If / ~ —1.0, the 
2 ++ state may also exists. From our numerical calculation, it mass should be about 

M = Mjy + M D , - E b ^ 4008.5MeV. 
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TABLE III: Binding energy Ej, (in MeV) of molecular states vs paramter /. 



/ 


-0.8 


-1.0 


-1.5 


-2.0 


-3.0 


(DD)0 ++ 






11.8 


28.1 


71.7 


{DD*)l ++ 


1.8 


6.2 


24.2 


48.7 


109.3 


(D*D*)2 ++ 




5.4 


21.4 


43.5 


98.1 


(DcDc)0 ++ 

\ J -'S- L -' s ) w 










5.3 


(D S D*)1++ 










15.4 


(D*D*)2 ++ 










14.0 


(BB)0 ++ 


8.8 


15.3 


35.2 


58.5 


110.9 


(BB*)l ++ 


16.8 


26.8 


55.8 


88.5 


160.4 


(B*B*)2 ++ 


16.0 


25.5 


53.4 


84.7 


153.6 


(B s B s )0++ 






3.4 


9.3 


25.7 


(B S B* S )1 ++ 


- 


1.9 


8.9 


18.9 


44.4 


(B* S B* S )2++ 




1.8 


8.5 


18.2 


42.9 




TABLE IV: rms 


(r 2 ) 1 / 2 (in fm) of molecular states vs paramter /. 




f 


-0.8 


-1.0 


-1.5 


-2.0 


-3.0 


(DD)0 ++ 






1.33 


1.00 


0.76 


(DD*)1++ 


2.67 


1.68 


1.07 


0.86 


0.69 


(D*D*)2 ++ 




1.79 


1.12 


0.91 


0.72 












1.63 


(D S D*) 1++ 










1.12 


(D*D*)2++ 










1.17 


(55)0++ 


1.11 


0.94 


0.74 


0.64 


0.54 




0.92 


0.80 


0.66 


0.58 


0.50 




0.94 


0.82 


0.67 


0.59 


0.51 


(B s B s )0++ 






1.42 


1.01 


0.74 


(B S B* S )1++ 




1.77 


1.03 


0.81 


0.63 


(B* S B* S )2++ 




1.81 


1.05 


0.83 


0.64 
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The Belle collaboration has reported a X(3930) state of 2 ++ at mass M = 3929 MeV 26], 
which is a candidate of the cc charmonium 2 3 P 2 excited state x' C 2- However its mass is 
aboult 40 MeV below the quark model calculations (in the Bhaduri model, the x' C 2 m ass is 
M = 3963.53MeV) . Several authors have discussed the mass shifts of the coupled channel 



effect 



27 



28 



29 



301 ] . If the above 2 ++ molecular state exists, the final state interaction will 



be important. 

Now, we turn to the possible molecular states in other sectors. We observe that in 
DgDg sector, there are no such molecular states due to the mass increase of m s of the 
light quark pair. On the other hand, in the _£?(*)_£?(*) sector, the binding energy increased as 
the mass increase of rrib of the heavy quark pair. So these bottomium-like molecular states 
with 1 ++ , 2 ++ , and even ++ may also exists. We can also observe the Bg Bg molecular 
states 1++ and 2++ if / 1.0. 



IV. SUMMARY 



The quark model is extended by introducing the s-channel one gluon exchange interaction. 
The interaction has no effect on the inner quark structure of conventional qq mesons and 
qqq baryons. Since the interaction is short ranged — a smeared (^-interaction, the effect on 
the long ranged hadron-hadron interaction is expected to be very small. So the significant 
effect of this interaction is only in the so called hidden flavor states of multi-quark system. 

We have calculated the heavy quark molecular states of qQQq with Q = c,b and q = 
u, d, s. We find that the interaction can be strong enough to bind the 1 ++ and 2 ++ states, and 
possibly ++ states. To compare with the recent experiments, the X(3872) is a candidate of 
the DD* molecular state and the V(3930 may be the x' C 2 state which couples to the D*D* 
2 ++ molecular state. The calculation shows that it is more easy to bind the bottomium-like 
molecular states. Thus we expect that the similiar bottomium-like molecular states of 1 ++ 
and 2 ++ should also exist. 
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